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Abstract 

A completely rigorous first-principles calculation of the charged particle stopping power has recently 
been performed by Brown, Preston, and Singleton (BPS). This calculation is exact to leading and 
next-to-leading order in the plasma number density, including an exact treatment of two-body 
quantum scattering. The BPS calculation is therefore extremely accurate in the plasma regime 
realized during the ignition and burn of an inertial confinement fusion capsule. For deuterium- 
tritium fusion, the 3.5 MeV alpha particle range tends to be 20-30% longer than most models in the 
literature have predicted, and the energy deposition into the ions tends to be smaller. Preliminary 
numerical simulations indicate that this increases the pR required to achieve ignition. 
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I. INTRODUCTION 



The charged particle stopping power of a hot plasma plays a critical role in whether 
an inertial confinement fusion (ICF) capsule will undergo thermonuclear burn and ignition. 
For a robust experimental setup, in which the laser energy is sufficiently high to assure 
ignition and full burn, the fine details of the stopping power are not so relevant; however, in 
more marginal settings on the threshold of ignition, these details are likely to play a more 
important role. Conventional wisdom states that our current knowledge of the stopping 
power is probably "good enough"; however, I will present some analytic and numerical 
results to illustrate that the subleading order physics neglected or missed by most stopping 
power models, the terms of order unity inside the Coulomb logarithm, can in fact lead to a 
noticeable effect on ignition. 

In this study, I would like to explore some implications of the stopping power calculation 
of Brown, Preston, and Singleton (BPS) This calculation of the charged particle stopping 
power, which includes the energy splitting of the projectile as it slows down and a rigorous 
treatment of the quantum to classical transition, is near-exact for the weakly coupled plasmas 
relevant for ICF ignition. I should also emphasize that the BPS stopping power is not a model 
of the Coulomb energy exchange, but a controlled first principles calculation of this process. 
For a detailed pedagogical explanation of the technique, please see Refs. 0, S, 3- We may 
therefore study the effects of stopping power on ignition and burn in a systematic fashion, 
and unlike utilizing a model, we know the range of validity of the BPS stopping power, along 
with an estimate of the error in any given plasma regime. Rather than launching straightway 
into the simulation of real ICF capsules, with their concomitant computational and physical 
complexities, I shall instead start more systematically and investigate the simpler problem 
of the deuterium-tritium (DT) microsphere studied by Fraley et al. oj. We shall find that, 
for these so called Fraley spheres, the BPS stopping power increases the ignition threshold 
significantly compared to typical models of the stopping power found in the literature. For 
example, the value of pR necessary to attain a self-sustaining burn in a Fraley sphere increases 
by 10% in comparison to the well-known state-of-the-art model of Li and Petrasso [6|. 

Specializing to the case of DT fusion, the dominant boot-strap heating mechanism nec- 
essary to achieve a self-sustained burn is a particle energy deposition in the background 
plasma. All things being equal, the larger the stopping power dE/dx, the greater the energy 
deposition of the supra-thermal 3.5 MeV a particles, and therefore the easier it becomes to 
achieve a self-sustained burn through a particle heating of the ions. In other words, a larger 
stopping power leads to a more efficient self-heating, and conversely, a smaller charged par- 
ticle stopping power would tend to make ignition harder, since this would increase the range 
of a particles and lead to smaller energy deposition per unit volume. Other factors, such as 
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the relative amount of energy that the a particle apportions between plasma electrons and 
ions as it slows down, which is determined from dE e /dx for electrons and dE l /dx for ions, 
can also be quite important, with larger ion heating being more favorable for ignition. 

Let us start by reviewing a study performed by N.M. Hoffman and CL. Lee Q, in which 
they compared the stopping power model of 'Corman-Spitzer' (CS) [8] to that of C.K. Li 
and Petrasso [6]. A significant difference between these two models is that CS is valid only 
in a velocity window in which the projectile speed is much greater than the thermal ion 
velocity but much smaller than the thermal electron velocity, while Li and Petrasso (LIP) is 
more general and was constructed to remain valid for all non-relativistic velocities. The LIP 
stopping power tends to be smaller than that of CS, and consequently the a particle range 
for LIP is somewhat larger than for CS. By the aforementioned considerations, we would 
therefore expect the ignition threshold for LIP to be greater than that predicted from CS. 
Interestingly, however, Hoffman and Lee found that both models produce the same ignition 
threshold for the DT microspheres of Ref. [5], despite the longer range of the a particle 
predicted by LIP tJ. In addition to this, a preliminary study I have begun indicates that 
other models, such as that of Ref. j^] , produce an ignition threshold almost identical to that 
of CS and LIP. As pointed out by Ref. 0], these observations can be explained by noting 
that the models of CS and LIP differ most near the low energy thermal regime of the a 
particle projectile, while their stopping powers are almost identical at higher energies near 
the 3.5 MeV threshold (at even higher energies, CS and LIP begin to diverge, but this is 
well above the production threshold). Consequently, the stopping power models of CS and 
LIP act to slow the a particle down in almost exactly the same way throughout most of 
its history, since they give equal values for dE/dx throughout most of the projectile energy 
regime, and this is why the ignition profiles are virtually identical. In contrast, the stopping 
power of LIP and BPS differ significantly even in the high energy regime traversed by the 
a particle, and as we shall see, this produces a marked difference in the ignition profiles of 
LIP and BPS. 

Hoffman and Lee also observed that certain diagnostics, those sensitive to small projec- 
tile energies of order the thermal background, are quite different between CS and LIP. For 
example, Hoffman and Lee found LIP to be in agreement with experiment for the spectrum 
of fastprotons in a D 3 He filled capsule implosion, while CS disagreed with the experimental 
data 7, [lo| - The diagnostic signatures explored by Hoffman and Lee tend to be more sensi- 
tive to the low energy region of dE/dx, where CS and LIP have the largest discrepancy, and 
this explains the differences in diagnostics [7]. However, for diagnostics that are sensitive 
to the soft regime of the projectile, one must compute dE/dx extremely accurately. This 
presents a potential problem for models that are inaccurate near the thermal regime, since 
their diagnostic signatures are sensitive to this regime. However, as explained in Section II 
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of Ref. this criticism does not apply to the BPS stopping power, and consequently BPS 
should provide quite reliable diagnostics for such processes. 



II. CONTEXT 

The basic physics responsible for charged particle stopping in a plasma is energy ex- 
change through Coulomb interactions, a process that is quite similar to the Coulomb energy 
exchange that drives the temperature equilibration between plasma components at differing 
temperatures. Both of these processes, stopping power and temperature equilibration, are 
local, in that they do not involve energy transport in space. This should be contrasted with 
thermal conductivity, in which heat flow by electrons or ions is not only driven by Coulomb 
energy exchange, but also by the transport of energy from one spatial location to another. In 
this sense, the basic physics of the stopping power is identical to the physics of temperature 
equilibration between plasma components, both of which differ somewhat from the physics of 
thermal conductivity. I will have more to say about thermal conductivity in a future section, 
but for now, let us concentrate on the rate at which the charged particle loses energy as it 
traverses the plasma. For the Coulomb potential, a straightforward calculation of the energy 
exchange rate diverges logarithmically in both the short and long distance regimes [llj ]. To 
obtain a finite result, we must introduce ad hoc short and long distance cutoffs b min and 6 max 
as we integrate over the impact parameter b, in which case the rate takes the form 



The coefficient K is an exactly calculable analytic prefactor that depends critically on the 
physical process under consideration, while the logarithmic term is called the Coulomb loga- 
rithm. The exact values of the short and long distance cutoffs can only be estimated, and this 
is a dominant source of non-systematic error in model building. Choosing the values of b m \ n 
and 6 max based on physical considerations (however well motivated they might be), rather 
than direct calculations from theory, is what I will call model building. We may increase the 
sophistication of models such as Eq. (12.11) by including collective effects or improved short 
distance physics, as in Ref. [§], but I will still refer to this as modeling. In many cases in 
plasma physics, modeling is essentially the only means by which to proceed, particularly in 
complicated cases like warm dense matter. In some situations, however, such as the weakly 
coupled plasma of a burning ICF capsule or the hot plasma at the center of the sun, we can 
avail ourselves to simple perturbative analytic techniques rather than model building, and 
this is the path taken by BPS. 





Coulomb Logarithm (CL) 
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Before turning to the rigorous calculation of the stopping power, let us investigate further 
the model building methodology and the physical arguments used to choose the short and 
long distance cutoffs. Debye screening sets the scale for the long distance cutoff, and we 
expect on physical grounds that 6 max = c M k _1 with the dimensionless coefficient c M ~ 1, 
where k is the Debye wavenumber. The coefficient c M is by no means a constant, but rather, it 
is typically a function of the various plasma parameters (and the various projectile parameters 
in the case of stopping power), which I will write schematically as c M = c M (m,T,n). Short 
of a rigorous calculation, there is no prescription for determining the value or functional 
form of c M , and indeed, choosing the functional form of c M is part of the model construction 
(usually c M is set unity). There are also various choices for the Debye wavenumber k; for 
example, should one use the electron screening length k~ 1 or the total screening length k^ 1 
determined from the electrons and ions? The answer to this question is of course process 
dependent; for example, if the ions can be treated as static and screened by the electrons, 
then we should take 6 max = c M k^ 1 . However, to be rigorous, such choices must come out of 
the calculation rather than being put into the calculation by hand. The situation for 6 min 
is even less clear. The short distance cutoff is set by scattering. In the extreme classical 
regime, we take 6^ in = x (distance of closest approach), and in the extreme quantum 
regime we take 6^ in = x (momentum transfer)/^. In either case, the functions and 
are of order unity, but absent a rigorous calculation, this is all we can really say about them. 
To interpolate between the extreme quantum and classical regimes, one usually employs 
an ad hoc scheme such as 6 min = [(&min) 2 + (^min) 2 although in Ref. {^J it was shown 
that this simple interpolation is missing an important logarithmic contribution. From these 
considerations, it is clear that the model building process is at best logarithmically accurate, 
that is to say, that the coefficient inside the Coulomb logarithm is known only to within a 
factor of order unity, and a complete and rigorous calculation is the only way to really settle 
the issue (this is why it is not surprising that the coefficient inside the logarithm varies by 
an order of magnitude within the literature, and with few exceptions, there is no reason to 
prefer one model over another 1 ). 

Another point worth stressing is that the Coulomb logarithm itself is process dependent, 
despite the fact that the same basic physics of Coulomb energy exchange is at work in the 
stopping power and temperature equilibration. For example, if we consider the charged 

1 There are calculations in the literature, such as that of Gould and DeWitt [jjj], that do obtain the correct 
coefficients inside the logarithm. However, to my knowledge, none of these calculations are systematic, 
in that they retain spurious higher-order terms and do not provide an error estimate. Therefore, without 
independent verification such as BPS, these calculations do not have the power to know just how accurate 
they are. 
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particle stopping power and the electron-ion temperature equilibration rate, we can write 

^ = K p \n\ b ^\ (2.2) 
dx I o min I 

dS PJ 



dt 



} ( T <- T >)- 



K„lni-J^UT,-T, . (2.3) 



The first expression (12.21) represents the energy loss per unit distance of a projectile [the 
charged particle stopping power], while the second expression ( 12. 31) is the energy exchange 
rate per unit volume between electrons at temperature T e and ions at temperature T x [this 
rate is proportional to the temperature difference, which I have explicitly indicated in 
Eq. (I2.3p ]. I have placed a bar over the short and long distance cutoffs in Eq. (12.31) to 
indicate that the Coulomb logarithm need not be the same as for other processes. Con- 
centrating on the charged particle stopping power from here on, I will rewrite the generic 
stopping power model (I2.2p in the form 

^^model 

— = K cs lnA model . (2.4) 

dx 

Here, I am emphasizing that the leading order coefficient K cs is exactly known [sj], while the 
coefficient A mo d e i inside the logarithm is only known to an order of magnitude. The reason 
I do not call the BPS calculation a model is because Ref. [lj calculated the terms under the 
logarithm exactly from first principles, including the quantum to classical transition, along 
with a controlled estimate of the error, which, in contrast to Eq. (12. 4p . I will write as 

dE BPS 

— - — = K cs lnA BPS + y controlled jgrror term / . (2-5) 

small in weakly coupled plasma 

As I will discuss more fully in the next section, any thermodynamic quantity in a plasma 
can be expanded in integer powers of a dimensionless plasma coupling parameter g [131 ] . 
Measuring temperature in energy units and taking the electrostatic units to be rationalized 
Gaussian, the choices I will employ throughout this paper, the coupling takes the form 2 

e 2 n 1 

« = (2 - 6) 

where k is the Debye wavenumber. The ^-expansion admits nonanalytic terms such as \ng, 
and indeed, such terms are essential in capturing the interplay of short and long distance 



2 In ordinary nonrationalized Gaussian units we would write g = e 2 n/T . There is an independent coupling 
parameter for each plasma component, and this must be taken into account for a real calculation; however, 
for the purposes of explanation, we may consider an isolated plasma component. Finally, the usual plasma 
coupling r, defined in terms of the interparticle spacing, is related to the expansion parameter by g 2 oc T 3 . 
Integer expansions in g are therefore expansions in fractional powers of T. 
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physics. As we shall see in the next section, the stopping power may be systematically 
expanded in the form 

dE BPS 



ilr = -Ag^ + Bg^ + 0(g s ) , (2.7) 

LO NLO 

where I have indicated the leading order (LO) and the next-to-leading order (NLO) terms. 
To get a feel for the size of g in a plasma, at the center of the sun we find g = 0.04. 
In a weakly coupled plasma in which g <C 1, the error terms, which I have denoted by 
0(g 3 ), are quite small and the expansion (12.71) will be near-exact, provided of course that we 
know the coefficients A and B exactly. The coefficient A = A(m, T, n) is well known, while 
B = B(m, T, n) was calculated by BPS [1], and as a matter of completeness, the final results 
of the BPS calculation are displayed in Appendix|A] To make the connection with expression 
( 12. 5j) . I will write the leading order coefficient as K cs = Ag 2 , and define the dimensionless 
coefficient C = exp{— B/A}. We can then express the rate (12.71) in the form 

dE BPS 

— — = K cs \nA BPS + 0(g 3 ) , with lnA BPS = -\n{Cg} . (2.8) 

ax 

This gives the exact Coulomb logarithm since we can calculate C = C(m, T, n). Incidentally, 
and this cannot be overemphasized, this perturbative methodology obviates the need for the 
cutoff parameters b min and & max , and they can be dispensed with from here out. The long 
and short distant cutoffs are really only part of a heuristic device that, in my opinion, often 
introduces more confusion than it purports to settle. 

As discussed at length in Ref. to leading and next-to-leading order, and only to this 
order, can we decompose the stopping power into its contribution to the various plasma 
species, 

^ = £r§T = E. + ^> • (2.9) 

where we identify dE PPS /dx as the stopping power contribution from species b. The notion 
of dividing the energy into contributions uniquely associated with individual plasma species 
is valid only to order g 2 , as three-body and higher collective effects render this division 
meaningless when working to order g 3 and higher. For the BPS stopping power, the quan- 
tity dE PPS /dx depends not just upon parameters associated with the 6-species, but upon 
the plasma parameters of all other species through the dielectric function. This should be 
contrasted with the corresponding quantity dE\ lp fdx for LIP (except for the Coulomb log- 
arithm In A" p ), which depends only upon the plasma conditions of species b. Since we shall 
be comparing BPS with LIP, I will close this section by presenting the LIP stopping power 
in some detail. 
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Li and Petrasso [6j modeled dE/dx by combining a generalization of the Fokker-Planck 
equation to account for short-distance collisions and a well-chosen term involving a step 
function to include long-distance collective effects. They define a Coulomb logarithm by 
using a minimum classical impact parameter that interpolates between the classical and 
quantum regimes, as described earlier. Using rationalized cgs units for the electric charge, 
the LIP stopping power for species b can be written 



dE L 



dx 



dE^ IP 
b dx 



hA-TT 



hi, 



6 47r I3 h m b vl 



P L 



G ( - Pbiribvl In A 



(2.10) 



where the contribution to species b is dE b IP /dx. Here, the function multiplying the Coulomb 
logarithm is defined by 

r "" (2-11) 



where 



and 



G(y) 



m b d 
m p dy 



nip L 



1 d ' 
dy. 



dzz 1 ' 2 , 



H(y)+6(y-l) In (2e"V /2 ) 



(2.12) 



(2.13) 



with 6{x) being the unit step function: 6(x) = for x < and 6(x) = 1 for x > 0. Li and 
Petrasso define a Coulomb logarithm in terms of the combination of classical and quantum 
cutoffs as described above, namely 

1 



In AT 



where the short distance cutoff is defined by 



(2.14) 



(2.15) 



2m pb u b J ' \4imip b ul) 
in which m vb = m p mb/(m p + nib) is the reduced mass of the projectile and species b, and 

2 



u,, 



v 2 p + 



(3 b m b 



(2.16) 



defines an average of the squared projectile and thermal velocities. In rationalized units, the 
Debye wave number is 



h'i, 



el n b 



(2.17) 



In future sections, we will also need the plasma frequency, and in rationalized units it takes 
the form 



el n b 
m b 



(2.18) 
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III. ANALYTIC CONSIDERATIONS: THE BPS CALCULATION 



In this section I will review the salient features of the BPS calculation This calculation 
includes both hard (short distance) physics and dynamic collective (long distance) physics, 
joined together exactly and unambiguously, and systematized by a power series expansion in 
the plasma coupling constant g. Ref. [ij calculates both the charged particle stopping power 
and the electron-ion temperature equilibration rate, thereby providing an exact calculation 
of the coefficient inside the Coulomb logarithm for both processes in ( 12.2ft and ( 12.31) . An ad- 
ditional feature of the BPS calculation is that it also provides an exact interpolation between 
the extreme classical and extreme quantum regimes. The calculation exploits a procedure 
in quantum field theory known as dimensional regularization, or dimensional continuation 
as I will call it here. The basis idea is that divergent theories exhibit finite poles of the 
form l/(v — 3) when analyzed in an arbitrary number of spatial dimensions v. We can then 
manipulate finite quantities in such a way that we preserve the delicate relation between 
long and short distance physics. In a physical process, the divergent pole terms cancel and 
we can set the number of spatial dimensions to v — 3 at the end of the calculation, thereby 
giving a finite result with correct long and short distance physics. 

To start the discussion, let x and v denote the //-dimensional position and velocity vectors 
of a particle. The Coulomb potential for two particles separated a distance r = |x — x'| 
is V v (r) = C v e 2 /r u ~ 2 , where C v = Y[u/2 — 1)/Att u ^ 2 is a spatially dependent geometric 
factor. This potential follows directly from a simple multidimensional generalization of 
Gauss' Law. For every species b, the single-particle distribution function /& will be defined 
so that /b(x, v,t) d v xd v v gives the number of particles of species b in a small hypervolume 
d v x about x and d u v about v at time t. We shall take the plasma components to be 
Maxwell-Boltzmann, although with more work, the situation can be generalized to Fermi- 
Dirac statistics as well [3J. In the case of a projectile p, the distribution function f p will be 
peaked about a specific point in phase space, and the stopping power can then be calculated 



Expression (13.11) is, of course, extremely problematic in three dimensions, but completely 
finite when v ^ 3. We can gain some additional insight into the divergence problem, however, 
by returning to three dimensions for a moment. In three dimensions, if we knew the exact 
form of f p in the background of the other distribution functions /& of the various plasma 
components, then Eq. (13. ip would be completely finite and well defined in three dimensions. 
Unfortunately, the requisite solution for f p can only be obtained by solving the full BBGKY 



by 




(3.1) 
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hierarchy of kinetic equations, an impossible feat. The divergence problem in Eq. (13.11) 
arises only because we must use an approximation for f p , usually obtained by truncating 
the BBGKY equations to a first-order kinetic equation, such as the Boltzmann or Lenard- 
Balescu equation, and these approximates miss correlations that would otherwise render 
Eq. (13.11) finite in three spatial dimensions. Curiously, these divergences occur only for the 
Coulomb potential, and only then in three spatial dimensions! 

These observations provide a path forward. Let us return to an arbitrary number of di- 
mensions, where the Coulomb potential provides finite results. Let us also define multipoint 
correlation functions in a similar manner to the /'s, and in this way we can construct the 
full BBGKY hierarchy in an arbitrary number of dimensions. For simplicity I will drop the 
subscript on the distribution functions. In dimensions v > 3, the standard textbook deriva- 
tion of the Boltzmann Equation (BE) goes through without the infrared divergent scattering 
kernel found in three dimensions, and I will denote this finite i/-dimensional scattering ker- 
nel by the shorthand notation B v [f]. The important point to emphasize here is that the 
v- dimensional Coulomb potential V u ~ e 2 /r u ~ 2 emphasizes short-distance over long-distance 
physics when v > 3, and this means that the BBGKY hierarchy reduces to the Boltzmann 
equation to leading order in g in these dimensions: 

df 

BBGKY =^^ + vVJ= B u [f] to LO in g for v > 3 . (3.2) 

Here, the i/-dimensional spatial gradient has been denoted by V x . In dimensions v < 3, the 
standard textbook derivation of the Lenard-Balescu equation (LBE) goes through without 
the ultraviolet divergent scattering kernel found in three dimensions, and I will denote this 
finite i/-dimensional scattering kernel by the shorthand notation L v [f]. In dimensions v < 3, 
the Coulomb potential V v (r) emphasizes long-distance physics over short-distance effects, and 
consequently, to leading order in g, the BBGKY hierarchy reduces to the Lenard-Balescu 
equation (LBE) in this spacial regime: 

df 

BBGKY 4-^+v-VJ = L u [f] to LO in g for v < 3 . (3.3) 

Space does not permit us to write down the exact forms of B v \f\ and L v [f] here, but one 
may consult Ref. [l| for these expressions. These kinetic equations allow one to calculate 
the stopping power in v > 3 and v < 3, the results of which are presented in Sections 8 and 
7 of Ref. l|, respectively. These calculations involve performing a series of momentum and 
wave number integrals in arbitrary dimensions u, and they take form [l] 

dE > 1 r n 2 a 2 

V -— d v p ^— B v [f) = H(u) + 0(v - 3) : LO in g when v > 3 ,(3.4) 



dx v p J 2m p 

dx v r , 



[ d u p ^— L v [f] = G{v) ^— + 0(3 -v) : LO in g when v < 3 . (3.5) 
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The analytic expressions for H{y) and G(v) are rather complicated, 3 and space does not 
permit their reproduction here. In this paper, we are only interested in their analytic prop- 
erties as a function of v. In particular, the coefficients H{y) and G(v) can be expanded in 
powers of e = v — 3, and we find 

H{v) = -A + eH 1 + 0{e 2 ) and G(v) = -A + e G x + C(e 2 ) . (3.6) 

For our purposes, we do not require the exact forms of Hi, G\, nor that of the leading 
term A. It is sufficient to note that the leading terms in Eq. (13.61) are equal, so that 
H{y = 3) = G{y = 3). This is a fact that arises from the calculation itself, as it must, and 
it should be emphasized that this equality is not arbitrarily imposed by hand. It is a crucial 
point that the leading terms are identical, as this will allow the short- and long-distance 
poles to cancel, thereby giving a finite result. 

Since the rates dE^ /dx of Eq. (13.41) and dE^ /dx of Eq. (13.51) were calculated in mutually 
exclusive dimensional regimes, one might think that they cannot be compared. However, 
even though Eq. H 3 . 5 f) was originally calculated in v < 3 for integer values of u, we can 
analytically continue the quantity dE< / dx to values u > 3 (for simplicity I have been omitting 
the functional dependence on v from the stopping power, but it is there nonetheless). This 
process of analytic continuation to non-integer values of dimension is analogous to way in 
which the factorial function n\ on the positive integers can be generalized to the Gamma 
function T(z) over the complex plane, including both the positive and negative real axes. We 
can then directly compare Eqs. (13 .4p and (13. 5p . Upon writing the ^-dependence of Eq. (13.51) 
as g 2+ ( u ~ 3 \ when v > 3 we see that Eq. (I3.5P is indeed higher order in g than Eq. (13 .4p : 

dE< 2+(z,-3) 



V 



dt 



G(u) + 0(u - 3) : NLO in g when v > 3 . (3.7) 



By power counting arguments, no powers of g between g 2 and g v ~ x can occur in Eq. (13.41) 
for v > 3, and therefore Eq. (13.51) indeed provides the correct next-to- leading order term 
in g when the dimension is analytically continued to v > 3. The individual pole-terms in 
Eqs. (13 .4p and (13.71) will cancel giving a finite result when the leading and next-to-leading 
order terms are added. The resulting finite quantity will therefore be accurate to leading 
and next-to-leading order in g as the v — > 3 limit is taken: 



dE p 

— — = hm 

dx i^3+ 



dE> dE< 
+ 




+ 0(g 3 ) . (3i 



3 For the related process of electron-ion temperature equilibration, in contrast, the expressions for H {y) and 
G{v) are quite simple. 
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FIG. 1: A comparison of the stopping power of BPS to Li and Petrasso (LIP) for a 3.5 MeV a 
particle in a equimolar DT plasma. The plasma temperature is T = 3 keV and the electron number 
density is n e = 1 x 10 25 cnr 3 , giving a coupling g = 0.01. For BPS, the range and the energy 
deposited into ions and electrons are R* PS = 30/um, Ef ps = 0.38 MeV, and £* PS = 3.16 MeV, 
respectively. For LIP, these corresponding quantities are -R" p = 25fim, E^ IP = 0.45 MeV, and 
££ip = 3.09 MeV. The percent differences between BPS and LIP are: AR a = +20%, AE l = -16%, 
and AE e = +2%, where AX = (X BPS - X LIP )/X LIP . 



Note that this does not lead to any form of "double counting" since we are merely adding the 
next-to-leading order term (13.71) to the leading order term (13.41) at a common value of v > 3. 
We are now in a position to evaluate the limit in Eq. (I3.8p . Defining e = v — 3 as before, 
note that g e = expjelng} = 1 + elng + 0(e 2 ), which gives the relation 



9 - = - + \ng + 0{e) . 
e e 



(3.9) 



Substituting Eq. (13. 9ft into Eq. (13.71) . adding this result to Eq. (13.41) . and then taking the 
limit gives 



dE p 
dx 



-Ag 2 \ng + Bg 2 + 0(g ? - 



(3.10) 



with B = Hi — Gi, in agreement with Eq. (12. 7p . In this way, BPS has calculated the 
charged particle stopping power accurate to leading order and next-to-leading order in g. 
For completeness, the full BPS stopping power is presented in Appendix [A] 

As alluded to earlier, the BPS calculation predicts the range of the 3.5 MeV a particle in 
a hot DT plasma to be 20-30% longer than typical plasma models in the literature with a 
smaller energy deposition into the ions. For LIP and BPS, this is illustrated Figs. [T]and[2J 
As we shall see in the next section, the longer a particle range and less efficient ion heating 
tend to make ignition more difficult to achieve. 
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LP 

- BPS 




x (/Ltm) 

FIG. 2: As in Fig. [Q BPS vs. LIP for a 3.5 MeV a particle, in a DT plasma at T = 30keV and 
n e = 1 x 10 27 cm -3 , giving a plasma coupling is g = 0.003; i?^ PS = 3.7/im, £\ BPS = 1.51 MeV, and 
E* FS = 1.51 MeV; R% p = 2.9/xm, £ : LIP = 2.15 MeV, and £ e LIP = 1.36 MeV; A_R Q = 28%, AS, = 7%, 
and AE C = 11%. 

IV. NUMERICAL CONSIDERATIONS: EFFECTS OF STOPPING POWER ON 
IGNITION 

Since this is a preliminary study, which I would like to keep as clean and simple as possible, 
I will not model a real ICF capsule here. Instead, I would like to look at the essential features 
of the stopping power without the complications of additional processes like hydrodynamic 
instabilities and thermal conductivity. I have already mentioned the latter problem, but 
I would like to make a few more comments. In a real ICF capsule, consistency requires 
that we model or calculate the thermal conductivity just as accurately in g as the stopping 
power. However, calculating the heat flow is harder than calculating the stopping power: 
one must not only contend with the short and long distance physics of the instantaneous 
Coulomb interactions, but one must also consider the rate of energy transport from one 
spatial location to another. It is not clear to me whether the Coulomb logarithm of the 
thermal conductivity can be calculated with the BPS methodology, since one must invert a 
heat kernel in addition to performing the appropriate multidimensional integrals. One must 
also consider the possibility of a particles emerging from the hot spot into the surrounding 
colder and more strongly coupled plasma, a region where the BPS stopping power may not 
apply. Including these effects is of course essential, but at this stage, they are somewhat 
premature. 

For the purposes of this study, it suffices to concentrate only on the time period after 
maximum compression and minimum volume, near the onset of thermonuclear burn. Stop- 
ping power plays less of a role in the implosion process, which is set in motion by an assembly 
of high intensity lasers or pulsed-power diodes that generate a radiation field that ablates a 
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FIG. 3: The yield in MJ as a function of initial pR in g/cm 2 for a DT microsphere of Ref. [H]. 
The initial uniform temperature is taken to be To = 3keV and the initial mass of the sphere is 
mo = 100 /Ug. The upper axis gives the initial density p in g/cm 3 . The dashed curve is the yield 
profile for the stopping power of Li and Petrasso (LIP), while the solid curve corresponds to that 
of Brown, Preston, and Singleton (BPS). The points labeled by (a) and (b) correspond to full yield 
1% maximum yield, respectively. The points labeled by (c) and (c') correspond to 10% maximum 
yield, the first for LIP and the second for BPS. The BPS stopping power has a longer a particle 
range than LIP, and it delivers less a particle energy to the ions. This has the effect of increasing 
the ignition threshold for BPS relative to LIP. 



material shell (such as beryllium or plastic) encapsulating the DT gas. In fact, I will ignore 
the presence of a surrounding shell completely, and instead look at the idealized case of 
Fraley et al. [5]. These authors consider an initially static microsphere of compressed DT at 
uniform temperature To and density po- Because of the non-zero pressure in the microsphere, 
the outermost surface of the spherical assembly acquires an outward velocity (comparable 
to the speed of sound in the DT plasma), which produces an incoming rarefaction wave. 
The net effect is that the center of the sphere remains slightly hotter than the surface, and 
consequently thermonuclear burn (if it occurs) will start at the center of the sphere and burn 
outward. The a-particle from the fusion reaction D + T — > a + n will provide the bootstrap 
heating that will initiate ignition and maintain thermonuclear burn. 

For definiteness, we will take the initial temperature of these so called Fraley spheres to 
be T = 3 keV, with an initial mass m = 100 fig. We shall vary the initial density p and 
measure the final yield Y produced in the simulation. To find the final state of the system, we 
numerically integrate the Navier-Stokes equations using a two-dimensional Lagrangian finite- 
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difference scheme. Radiative energy exchange is included via a multi-frequency diffusion 
approximation. Other transport processes are included as well, and plasma electrons and 
ions are each assumed to be characterized by separate Maxwellian thermal distributions, and 
exchange energy at a rate characterized by Spitzer theory. 4 The results of these simulations 
are illustrated in Fig. [3X where we plot the yield in MJ as a function of initial pR in g/cm 2 . 
The dashed curve is the yield profile predicted by the LIP stopping power, while the solid 
curve corresponds to that of BPS. The points labeled by (a) and (b) correspond to full yield 
and 1% of full yield, respectively. When full burn is assured, as in case (a), note that the 
yield is relatively independent of the stopping power. The points labeled by (c) and (c') 
correspond to the 10% maximum yield threshold for LIP and BPS respectively. As we see, 
the effect of BPS relative to LIP is to push the 10% threshold from (pR) UP = 2.2 g/cm 2 to 
(pR) BPS = 2.4 g/cm 2 , an increase in pR of approximately 10%. 



V. ELECTRON-ION TEMPERATURE EQUILIBRATION 

The final question I shall address is the effect of electron-ion temperature equilibration 
on the ignition threshold. The rate at which electrons and ions come into equilibrium, 

*k = -C el (T e -T 1 ) , (5.1) 

was also calculated in Ref. [lj to the same level of rigor as the stopping power, and this is 
presented in Appendix [B] in full generality. In the high temperature limit that applies to 
ignition, the rate coefficient takes the simple form 

uj} o I m P , . . , 1 



C PT = — ft? i I — ~- In A RP o , with In A 



o _ e \/ n rri BPS > vv.luj.1 '""bps — n 

III V Z7T i P L 



h 2 



M S \ - - 1 



(5.2) 



where 7 = 0.57721 • • • is the Euler constant, K e and u e are the electron Debye wave number 
and plasma frequency, and = J2i u i * s sum °f the squares of the ion plasma frequencies. 5 
Comparing Eqs. (15.21) and (12. 3p . we again see that there is no need to reference the heuris- 
tic scales 6 m i n and 6 max - The BPS calculation predicts a smaller rate than typical model 
calculations in the literature. This is ostensibly advantageous for ignition, since the ion tem- 
perature is more weakly coupled to the electron temperature and can run away more easily; 



4 The BPS electron-ion coupling rate calculated in Ref. [l[ was also used, and little difference between this 
near-exact result and Spitzer. 

5 Equation (|5.2[) corresponds to Eqs. (3.61) and (12.12) of Ref. [1], where I have taken this opportunity to 
correct a small transcription error: when passing from Eq. (12.43) to Eq. (12.44) in Ref. [l|, a factor of 
1/2 was dropped. Restoring this factor of 1/2 changes the additive constant outside the logarithm from 
the —7 — 2 that appears in Eq. (12.12) of Ref. [l[ to the constant —7 — 1 in Eq. (|5.2[) of this paper. 
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however, initial simulations indicate that the electron-ion equilibration has little effect upon 
ICF yield. Nonetheless, I would expect this rate to have some effect upon burn diagnostics, 
and I will explore this possibility in the future. 
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APPENDIX A: THE BPS CHARGED PARTICLE STOPPING POWER 

As a matter of completeness, I will present the full form of the BPS stopping power 
calculated in Ref. [l|. Consider a projectile of mass m p and charge e p moving through a multi 
component plasma at speed v p . Each plasma component is assumed to be in equilibrium 
with itself but not necessarily with the other components, and will be labeled by an index 
b (including the electron component). The mass of a component is m&, the number density 
is rib, an d the temperature is T&. Temperature will be measured in energy units, and I will 
use the notation /3b = for the inverse temperature. To second order in the plasma 
coupling g, the order to which we are working, the stopping power dE/dx of the projectile 
can meaningfully be broken into the contributions dEb/dx from the separate components, so 
that dE/dx = ^2 b dEb/dx. Each contribution can be expressed as a sum of three terms, 

dEb = fd_E^ + dE^\ + dEl 
dx \ dx dx ) dx ' 

where the first two arise from classical short and long distance physics, and the latter term 
from short distance two-body quantum diffraction. The stopping power dEb/dx is a function 
of the parameters of the projectile and the plasma background. For convenience, however, 
we shall emphasize the functional dependence only upon the projectile energy E p = ^m p v p . 
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The terms in Eq. (lAlj) are given by Eqs. (3.4), (3.3) and (3.19), respectively, in BPS 
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In the long distance contribution ( 1A3j) . we define # as the angle between the vectors k and 
v p , while F is equivalent to the plasma dielectric function, 



k 2 e(k,u = k ■ v p ) = k 2 + F(v p cos9) in which F(v) 
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In Eq. f 1A4|) the integration variable is 
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and the dimensionless quantum parameter is 
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The Diagamma function is defined by if) = T~ 1 dT /dz, so that 
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where 7 = 0.5572 • • • is the Euler constant. The total mass and the reduced mass are 

M pb = m p + m b — = — + —. (A10) 

m pb m p m b 

The sum of terms flA3[) and (IA2[) form the classical contribution, and the factor K is an 
arbitrary wave number that cancels in the sum of Eqs. (1A3I) and (IA2I) . It is convenient to 
set A" = Ac e . 



APPENDIX B: THE BPS TEMPERATURE EQUILIBRATION RATE 



For the multi component plasma described in Appendix [A] we assumed that the various 
components were in thermal equilibrium with themselves but not with each other. In prac- 
tice, however, components exchange energy via Coulomb interactions, and will equilibrate 
according to a common temperature according to 

d£ a b 



dt 



Cab (T a — T b ) . 



(Bl) 



As with the stopping power, the rate coefficient can be written as a sum of three terms, 

C-ab — {C a b,s + ^ab,lij + ^ab ' 

given by (12.31), (12.25), and (12.50) respectively in BPS Q: 
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The function F{v) in Eq. ( 1B4I) is defined by Eqs. ( 1A5I) and ( 1A6I) . and the strength of the 
quantum effects associated with the scattering of two plasma species a and b is characterized 
by the dimensionless parameter 
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where the square of the thermal velocity in this expression is defined by 

m a m b 
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In the limit (3 e m e <C /^m^ we have 



V 2 = — 



Upon taking the sum over ions, C el = J2fiei, and using the inequality f3 e m e -C /^m*, we can 
express the rate coefficient as [3] 

2 / a \ l/ 2 
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In the extreme quantum limit the term Aj in Eq. (1B10|) vanishes, and we obtain Eq. (15.21) . 
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